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ABSTRACT 


Calculations are made of temperatures vs depth and time 
in snow-free winter ice with hourly atmospheric heat fluxes as input 
data for the month of April 1968. Heat fluxes were derived from weather 
observations made at stations located in the Canadian Arctic and are 
assumed to be valid for the ice cover in western Parry Channel. A 
Simplified analytical model is justified on the basis that the prime interest 
is in temperature fluctuations of a daily period near the ice-air boundary. 
Results are shown for different assumed ice qinface albedos and short wave 


extinction coefficients. 
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INTRODUCTION 

The need to predict short-term temperature changes vs time in 
a sea ice layer arose out of the observation that surface tension cracks 
produced underwater noise and that the noise responded to decreases in the 
air temperature (Ganton and Milne, 1965). Of particular interest were the 
variations in ice temperatures vs time that existed during the period 
16 April to 27 April, 1968. During this period fluctuations in the air 
temperature corresponded to antiphase fluctuations in the underwater noise 
envelope recorded in M'Clure Strait. These fluctuations are shown in 
figure 1 and M'Clure Strait is shown in the map of figure 2. 

The calculation method used was to assume snow free winter ice 
of a given thickness and to synthesize the smooth interchange of heat 
fluxes with time between the ice and the atmosphere by the use of hourly 
step function changes in heating or cooling. By summing solutions of the 
equation of heat conduction appropriate to these step function changes 
in heat flux, temperaturesvs depth and time were obtained. 

Maykut and Untersteiner (1969) perform much more detailed 
calculations of a similar nature which were designed to predict the 
equilibrium conditions of the growth and decay of sea ice. Their input 
data was primarily mean monthly heat fluxes. Thus, given initial conditions 
of ice and snow thickness, the effect of various parameters could be 
explored with respect to the equilibrium conditions achieved after the passage 
of several years. In contrast, the calculations in this report use hourly 
heat fluxes as input data, the intent being to explore the maximum 
temperature changes (snow free) vs time that are likely in a sea ice layer 
with respect to variations in parameters such as the albedo and extinction 
coeeficient. Hourly heat fluxes, used for input data, were derived from 


nearby Arctic weather station observations - primarily from Mould Bay, 


AD ae 
for the period of interest (figure 2). 

The prediction of transient temperature changes in a floating 
sheet of sea ice is a prerequisite in the estimation of variations in 
lateral stresses and strains with depth. The compressive and tensile 
stresses which result from temperature changes are important in the con- 
sideration of ice pressures exerted on structures and in the occurrence of 
thermal fissures or cracks such as those described by Zubov (1943) and 
Kingery et al.(1963). Thermal tension cracks are important sources of 
impulsive underwater noises. Milne and Ganton (1971) summarize the charac- 
teristics of this type of under ice noise. 

It is difficult to predict the actual lateral stresses and 
strains induced by temperature changes in an ice sheet because of the 
general lack of knowledge of the macro-viscoelastic properties of sea ice. 
Tabata (1948,1967) shows that sea ice which is warm and saline is more 
plastic than when cold and less saline, and Voitkovsky (1967), in his 
measurements of the relaxation of stresses in fresh ice, shows that ice 
is a non-linear viscoelastic solid. The result is that sea ice has 
viscoelastic properties which depend on its temperature, salinity and 
stress history, all of which are unknown in detail. To complicate matters 
further measurements of salinities vs depth in winter ice of apparently 
uniform thickness vary greatly depending on the location as shown by 
Bennington (1967). 

In view of the above uncertainties it seems reasonable to make 
some simplifying assumptions in the calculation of temperatures vs time 
in a layer of sea ice. The most basic of these is to assume that the 


thermal conductivity and diffusivity are constants, independent of salinity 


ates 
and temperature. Other simplifying assumptions include a constant albedo 
at the ice-air interface and a constant extinction coefficient for the solar 
radiation which penetrates the ice and causes internal heating. Constant 
values for the conductivity , diffusivity and the extinction coefficient 
imply that the ice layer is thermally homogeneous so that the one-dimensional 
diffusion equation can be used, with boundary conditions that need only 
be applied at the ice-air and the ice-water interface. 
THE HEAT CONDUCTION PROBLEM 

Solutions of the equation of heat conduction 

Figure 3 shows the geometry of the one-dimensional heat conduction 

problem, where z = H is the ice thickness and z = 0 defines the ice-air 
boundary. The net heat flux, positive in the +z direction, is shown 


oe It is assumed that the net heat 


as QS, and is in units of cal cm_ 
flux, QS, enters the ice layer at z = 0, some of which heats the boundary 
layer at z = 0, while a radiative component is dissipated throughout the ice 
layer and the water beneath. The water is assumed to be a constant 
temperature bath. 

From p.29 of Carslaw and Jaeger (1960), the appropriate one- 
dimensional equation of heat conduction is: 
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where T = T(z,t) is the temperature at depth z and time t, 
k is the thermal conductivity, 
aoe a is the thermal diffusivity, with 
o the density of sea ice, 


c the specific heat; 


and A is the rate of heat production per unit depth. 


ay 
Since a net heat flux, QS, is assumed to have entered the ice 
at the plane z = 0, the boundary conditions are: 
(1) no heat conduction exists at z = 0 @TAz = 0 at z = 0), 
and (2ytaerz 7=°HST =P0. 
In a homogeneous solid subjected to radiant heating, the 
rate of heat production per unit depth decays exponentially so that: 


we A exp (-Kz) (2) 


where K is an “extinction coefficient". {if a part@oft the total neacterlux, 


QS, which causes internal heating,is termed, Q, then 


Q= At. exp (-Kz)dz = AU/k. (3) 

(Later, it is expedient to find a solution for heating which 
is concentrated near z = 0 by letting the extinction coefficient, K, 
approach infinity). 

Using (2) and (3), equation (1) becomes: 
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The total heat flux, QS, normally has several components, all 
of which can be assigned different extinction coefficients so that if 


QS = Q,*Q,+ ---Q, then (4) can be written as the set of equations: 
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where T= T)+T,+-.-T.. 

If all of one component of heat flux, say Qy is totally 
absorbed near z= 0, then T in (5), with the extinction coefficient K 
set to infinity, is also a solution. 

Considering the current problem, the total heat flux, QS, 
can be divided into two parts, one part QQ which produces internal heating 
and a second part Q, which produces surface heating at z = 0. 

The required solution to (4) for the step function of heat 
flux Q(t) = Q> t > 0 only, was obtained using the method on p. 130 of 
Carslaw and Jaeger (1960), where the temperature Tat bina) esez acall 
is expressed as the sum of a steady state and a time dependent solution of 


(4) as follows: 


T(z,t) = u(z) + w(z,t), (6) 
where 
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= ~ KQ, exp (-Kz)/k, u = u(z) only. 
Hence u(z) = - Q exp (-Kz)/ (kK) + Cz + C,> Cc, and C., constants. 


Applying boundary conditions <4 = 0 atez =20 “andju =, 0 atez =H; 
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u(z) = (Q,/k) 4 (H-z) - exp (-kH) exp (-K(H-z))-1] /K) (7) 


From (6) and (7), equation (4) now becomes: 
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a series solution to (8) is: 
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with w(z,0) = toe bcos ((2n41)1(2H)) | ; (9b) 
Since no heat is applied until t > 0, then from (6) 
u(z) + w(z,o) = 0, 

and from (7), 
w(z,0) = (Q,/k) | (-2)-exp(-Ki) [exp (-K(H-2)) -1] 7x | (10) 


The coefficients, ba in the sum (9b) are given by: 
H 
bo =11(2/H) Lo w(z,0) cos ((2n+1)1z/ (2H) )dz. (1) 


Substituting (10) in (11), and integrating: 
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From (6), (7), (9a) and (12), the solution to (4) for a step function 
of heating flux Q(t) = QQ t > 0 only, is termed T,(,t), where the 


subscript "I refers to internalsheating = andeiseetvenmpy : 
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x cos [ (2n+1)1z/ (2H) | | (13) 


For the case where all the heat flux causes heating solely 
in the vicinity of z = 0; a second solution, termed Ty(z,t), where the 
subscript "U" refers to upper surface heating, is obtained from (13) 


by letting the extinction coefficient, K, approach infinity; hence: 
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Figures (4a) and (4b) are plots of T, (@,t) from (13) and figure 
(4c) shows plots of Ty (Z5t) from (14). The curves in each case are for 
constant times, t, and the families are for 20 hr increments in time 
beginning at t = 0.5 hrs. In figure (4a), Qrur clive. £21 ace eand 
re Sequin aioe gt figure (4b), Q, = 1.75 cal ee Hore a Datet hepextinct ion 
coefficient K is decreased to 0.0085 cm. Figure (4c) shows the effect 
of K equal to infinity (surface heating only) and where Qy = 0.875 cal 
emer the ocher constants used in (13) and (14) were as follows: 


ice thickness, H = 200 cms; conductivity k = 17.5 cal ones Chr a 


and the diffusivity a~ = 30.2 cm*hr -. 
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Heat Fluxes 
Walmsley (1966) describes in detail the ice cover and surface 
heat fluxes in Baffin Bay and in particular the interelationship of surface 
weather observations with computed heat fluxes. His choice of expressions, 
describing the convective and latent heat fluxes,is used herein. 
| The heat flux balance can be expressed as follows: 
QR + QE + QC = QS, (15) 
where QR is the net radiation flux, 
QE is the latent heat flux brought about by the evaporation or 
condensation of water substance, 
QC is the sensible or convective heat flux, 
and QS is the net heat flux change stored in the sea ice 


and the ocean beneath. 


eS ae 
The convention adopted is that all heat fluxes are positive 
in the +z direction (see figure 3), i.e., when heat is added to the ice. 
Convective heat flux, QC 
Walmsley (1966) uses Shuleikin's (1953) formula for QC, where: 
Qc = 1.27 (T,-T) ‘cal em “hr ?, T eed (16a) 


or QC 0175 (T,-T)V cal cut hr? erp ae (16b) 


in (16) qT, is the air temperature exon) at a height of 8m, T is the surface 
temperature Co of the ice and V is the wind speed at 8 m height (m see) 
Latent heat flux, QE 
Sverdrup's (1951) formula is used by Walmsley, where the 


evaporation of water from a sea surface is given by: 


-l 
E = K (W.-W )V cn day) *, 


where Ws 5W, are water vapour pressures at the ice surface and in the air, 
V is the wind speed and 
Ku iSea coeftiilcient: 
The latent heat flux is given by 
QE = aLE 
where a is the density of the ice, 
L is the sum of the latent heats of fusion and sublimation and 
E is the evaporation. 
The water vapour pressures W; and Ne are expressed in terms of 
the ice surface temperature, Rac) and the dew point TC) at a height of 8m, 
and after establishing suitable values for the coefficient x, Walmsley 


uses the following expressions: 


Wee . [ exp (.0728T ,)-exp(.082T) ] jecalgcmpcin iG (17a) 
Foraves 6.2msec?; 
1.4 [exp(.0728T,)-exp(.082T)] v cal cm hr}, (17b) 


for V < 6.2msec!, 
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It can be seen that the heat fluxes QC in (16) and QE in (17) 
are sensitive to the ice surface temperature which is to be calculated. 
This difficulty is partly resolved by assigning the ice surface temperature 
calculated for the previous hour to the computation of QC and QE for the 
current hour. Both fluxes QC and QE interact at the ice-air boundary 
(z = 0 in figure 3). 

Radiative heat fluxes, QR 
Following the notation of Schwerdtfeger and Pounder (1963), 


the net heat flux from electromagnetic radiation is given by: 


where A. and AG are the total incident and outgoing radiative fluxes for 
all wavelengths. A. and Ao are separated into the short wave (0.2 to 4 micron) 
fluxes S. and So and the long wave (3 to 80 micron) fluxes Ls and Loe sO 
that A. = Si +L; and A, = Sot: The long wave radiation is entirely 
absorbed and emitted in the vicinity of the ice-air boundary (z = 0), 
so that Ly can be expressed as Ly = ele where ep is the long wave 
emissivity, o is the Stefan-Bobtzmann constant and Ty is the ice surface 
temperature (°K). The short wave radiation, on the other hand is partly 
reflected from and partly transmitted through the ice surface according to 
its albedo, a. Therefore Sy = as. and the short wave net heat flux is 
(1-a)S,. A fraction of (1-a)S, causes internal heating which changes with 
depth according to an "average" extinction coefficient, K, for the short 
wave radiation band. 

Maykut and Untersteiner (1969) introduce an internal heating 
flux I); where the difference between (1-a)S; and I, is a contributor to 


surface heating. It is difficult to estimate the fraction I¥/(Q-a)S; ) 


since it is likely to range from near unity for a low albedo to .1 or so 


a} Ope 
when the ice surface has a high albedo. This is equivalent to saying that 
there can be a thin layer near z = 0 which has a high extinction coefficient 
and at the same time reflects most of the short wave radiation. For the 
purposes of the calculations in this report, I, = -5(1-a)S,. To sum up, 

QR = L; - Lo - (1-a)S, (18) 
where the flux L. - L, + (1-a)S; - I) causes surface heating, and the flux I, 


causes internal heating. 


The Storage heat flux, QS 
QS is the sum of all three heat fluxes QR, QC and QE and appears 


as heat loss or gain in the ice and the ocean beneath. A small fraction of 
the short wave radiation is lost to the water, which is assumed to be an 
infinite energy sink. The balance is stored inz(l1) heat changes in the ice 
layer, (2) changes in the latent heat of fusion at the ice-water interface, 
and (3) in the convection of heat from the water. Lewis(1967) indicates that 
10% of this balance, under conditions of ice growth, comes from the water, 
while under equilibrium conditions, the remainder is devoted to the accretion 
of ice so that the latent heat of fusion is the primary sink of heat. 

For the purposes of the problem described in this report, the 
fraction of the short wave radiation energy lost to the water is ignored 
and it is assumed that 90% of the daily sum of QS, if negative, is devoted 
to the liberation of the latent heat of fusion in the accretion of ice. 
For each day, this accretion, in cms, is added to the initial thickness. 


The daily accretion is therefore assumed to be 


AH = - = QS en day"! (19) 


where <QS> is the daily average of (QC+QE+QR) in units of cal cm -“hr7?, 
The value for the latent heat of fusion, L, and the density of sea ice,p, 


at the growing interface were taken from Lewis (1967), where 


o aige 
=a) -3 : 
L = 67.5 cal gm = and op = .915 gm cm. Using these values, (19) becomes: 


AH = - .35 <QS>cm day. 


Heat flux data inputs 
Since no heat flux measurements were made over the ice cover 
of interest, reliance was placed on the weather observations made primarily 
at Mould Bay and at Resolute Bay. The intention was not to replicate true 
temperatures vs time in the ice sheet but to obtain typical temperatures 
which would be produced by similar weather systems during the month of April. 


Air temperatures, TT. dew points T,, wind speeds, V, and their directions 


d? 
were available from the bi-yearly "Arctic Summary" publications. for Mould 
Bay. These values are listed every 4 hours. For the purposes of calculation, 
hourly values were obtained by a linear interpolation procedure. Radiation 
observations were obtained from the ''Monthly Radiation Summary" publications. 
These are listed hourly, however not all necessary observations were 
obtainable from Mould Bay so that some of the missing radiation observations 
were obtained from Resolute Bay. 

The calculation of QC from (16) and QE from (17) is therefore 
relatively straight forward; however the calculation of QR in (18) from the 
available radiation information is somewhat more complex. The available 
radiation data of use was: 

RFl, the global solar radiation, downward and diffuse, 

RF3, the reflected solar radiation, upward and reflected, 

RF4, the net radiation that is the net flux of downward and 

upward total solar, terrestrial and atmospheric radiation. 


Therefore, at the weather station, 


4 
S,-S+(L,-L,) = RF1-RF3+(L;-e,0T,) = RF4, 


from which (L; -L,) = RF4+RF3-RF1. | (20) 
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‘i Available from the Department of Transport, Meteorological Branch, Canada. 
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than to assume that surface temperature T at the weather station approximated 


Recalling that Ly ae aT? = ea (T+273)", there was no other choice available 


T at the ice-air interface insofar as the long wave radiation was concerned. 
For the short wave radiation, RFl = S.5 so that over the ice 

S.-S) = (1-a)S; = (l-a)RF1, (1) 
where a is the albedo of the sea ice. Therefore from (20) and (21), 


(18) can be written as: 


QR = RF4 + RF3 - a(RFl) cal emschua.s (22) 
of which 

RF4 + RF3 - a(RF1) - 1, (23) 
causes surface heating, and 

I, = 0.5(1-a)RF1 (24) 


causes internal heating. I) is hereafter given the Fortran name 'RADQ'. 

The net heat flux which caused surface heating in the vicinity 
of z = 0 is given the Fortran name 'SURFQ' and is the sum of QC, QE and the 
radiation surface heating from (23): 

SURFQ = QR + QC + QE - RADQ. (25) 

As mentioned previously, values of QC and QE depended on the 
calculated surface-temperatures of the sea ice and were made available 
only for the hour that had passed. The radiation terms, air temperatures 
and dew points were stored in advance. 

Numerical procedure 

At any hour, the temperature vs depth depended on the past 
history of heat flux inputs, and since the solutions T, (,t) nel So) ano 
To (z,t) in (14) are applicable to step function heat flux inputs, then the 


temperature T at hour N was a result of a set of heat flux steps: 


aria 
(RADQ (N-m) - RADQ (N-m-1)), from (24) and 
(SURFQ (N-m) - SURFQ (N-m-1)), from (25), 
where mM SOOM 2ecieNe2 3 


The temperature vs depth at hour N can be expressed as follows: 


N-2 


1 ul 
TAZANJ evi (2,194 r1=0 T(z ,m+ 5) + T, (z,m+ 5) yz) 


with N = 2,3...NMAX, 

Qy in (14) set equal to equal to SURFQ(N-m) - SURFQ(N-m-1), 
and QQ in (13) set equal to equal to RADQ(N-m) - RADQ(N-m-1). 

T(z,1) is a linear temperature vs depth profile assumed as an 
initial condition. In (26) T(z,l) = (-1.56+H <QS(1)>/k) °C, where 
<€ QS(1) > is the daily average of QC+QE+QR in cal cm hr} for the first 
day and the sea water temperature is -1.56 ete 

Constants used in the calculations 

The thermal conductivity, k, and the diffusivity, ae = k/ (pc) 
are variables in sea ice which depend on the temperature and salinity 
variation with depth. Untersteiner's (1961) formulae express this 


variability as follows: 


k = 17.5 + 1.01 S(z)/T(z) cal em) °C7) pra), 


and 

Ao te k/ (pe) =k) (0N45240 aS 2d i (20) em hia, 
where S(z) is the salinity in parts per thousand (%) and T(z) is the 
temperature in ech 

Figure 5 b shows k and aa vs depth for the winter ice conditions 
shown in figure 5a and figure 5d shows k and Be vs depth for the polar ice 
conditions shown in figure 5c. The diffusivity is considerably more variable 
with temperature, as a result of the changes in the specific heat c. The higher 


specific heat is a result of the storage of heat in the changing volumes of 


214 
brine, which increases in volume nearer the melting point. A good 
description of this phenomenon is given by Schwerdtfeger (1962). 

The consequences of assuming that the conductivity, k, is 
a constant are relatively insignificant. The consequences of assuming a 
constant diffusivity, aan need, however, to be assessed in the light of the 
problem under consideration. For the steady state, ae: does not influence 
the temperature profile. When the heat fluxes vary with time, it can be 
seen from an examination of equations (13) and (14) that decreasing values 
of ae require increasing times t to obtain the identical temperatures, 
and therefore reflect the increasing heat energy storage capability of 
warmer sea ice. Consequently, an assumption of a constant value of att 
say 30.2, to represent reality for the upper 100 cm of winter ice (figure 5b) 
will predict transient temperature fluctuations which are increasingly 
too large as the ice-water boundary is approached. Because short period 
temperature fluctuations are rapidly damped with increasing depth, errors 
will only be significant with regard to longer period temperature fluc- 
tuations. Even these errors are likely to be irrelevent to the calculation 
of lateral tensile or compressive stresses because of the increased 
plasticity of sea ice with increases in salinity and temperature. 


For the purposes of the current problem, k = 17.5 cal en tPagdpyy-t 


and a* = 30.2, cm-hr 7}, 
Two other constants required are the albedo, a, and the 
extinction coefficient, K ciiaine The albedo is the ratio of the reflected 


to incident radiation and the extinction coefficient defines the decay 


in the intensity of the solar radiation flux with depth. Both a and K 


~m1Si- 

depend on the electromagnetic wavelength. Measurements of the albedo, 

a, are described by Langleben (1969), Briazgin (1959) and Maykut and 
Untersteiner (1969). From these authors, a value of a = .65 was chosen 

to be representative of the albedo of snow free winter ice in April. 

This value is likely to be, on the average, somewhat low, however the error 
is in the direction which would lead to higher surface temperature fluc- 


tuations. The extinction coefficient for the solar radiation penetration 


of winter ice does not appear to be specifically available in the literature. 


For sea ice, unspecified.,values of K are given by Thomas (1963) and 
Roulet (1970) where the latter author shows the variability in K with depth 
and electromagnetic wavelength. A value of K = .0170 was chosen to be 
representative of winter ice in April. 
RESULTS 

Figures 6a, 6b and 6c are results of calculations for the 
period 2300 hrs (LMT) 15 April to 2300 hrs (LMT) 27 April, 1968, inclusive. 
(Calculations actually commenced at 2300 hrs (LMT) on the 11th of April 
in order to diminish the effect of an incorrect thermal history). Figure 
6a shows predictions of ice temperature Elsi VS time act depths of 0:0; 
8.0 and 16.0 cms. Features shown are (1) a steady increase in the average 
temperature, (2) decreased amplitude and a phase shift with increased 
depth for the temperature fluctuations with a daily period, and (3) the 
low pass filtering effect of the ice on the higher frequency temperature 


fluctuations. Constants chosen were: conductivity k = 17.5 cal cm! clay 


diffusivity neh baci cmhr?, extinction coefficient K = 0.017 m7}, 
albedo a = 0.65. As in equation (24), one-half the short wave radiation 
flux was assumed to be absorbed at the ice-air boundary and the remainder 


internally. The initial thickness of the ice H was assumed to be 213 cm 


on the 12th of April. 


-1 
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Figure 6b shows the heat fluxes QR,QE and QC from (16). Glee 
and (22) in cal cm’“hr +. It can be seen that the latent heat flux, QE, 
is very small and that the radiative heat flux, QR, was dominant during the 
period 16 to 27 April, 1968. The en values of QR were near 1400 hrs, 
local mean time, where the ice cover heat gain was a maximum. The 
convective heat flux, QC, remained negative, indicating that the ice surface 
temperature remained warmer than the air temperature according to (16a). 
In fact, QC remained roughly in phase with QR indicating that radiation 
cooling from the bare ice surface at night resulted in ice surface tempera- 
ture decreases less in magnitude than air temperature decreases. Whether 
or not such a situation was realistic could only have been revealed if air 
temperatures over the ice had been available, rather than the air temperatures 
at Mould Bay. 

Figure 6c shows the internal heating flux, RADQ, from (24) 
and the surface heating flux, SURFQ, from (25) which were used in (26) 
to compute the temperature vs time in figure 6a. Daily variations in RADQ 
are a result of the presence and absence of solar radiation entirely. 

Calculations shown in Figures 7 and 8 are similar to figure 6 


except for changes in the albedo, a, and extinction coefficient, K, all 


other constants remaining the same. In figure 6 the calculations were for 


1 
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ah2,0265. and K,=40,017,cm,'; while in,fipure 7 for a < 0,75 and /K = 0/054 ene 
and in figure 8 for a = 0.55 and K = 0.0085 cmt, It was assumed that if 

the albedo was lower, the ice was likely to be more transparent, so that 

a lower extinction coefficient should be used; conversely, a higher albedo 
was assumed to be accompanied by a higher extinction coefficient. 


It is interesting to observe that the calculated amplitude of 


the daily fluctuations in ice temperature varied little between figures 7a 


ee f 

and 8a; although for the higher albedo in figure 7, the slowly rising trend 
in surface temperature was about 3°C lower as a result of more reflection of 
solar radiation. The daily periodic fluctuations in ice temperatures in 
figure 7 for depths of 0.0, 8.0 and 16.0 cms tended to remain similar to 


those in figure 8 as a result of the higher extinction coefficient used. 
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Figure 4a. T, (z,t) from equation (13), each curve being for a constant time, t, 


where t changes in increments of 20 hrs beginning at 0.5 hrs. 


Computations are for K-= 0.017 cm! and Q; = 1.75 cal cem~“hr) 
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Figure 4b. T, (z,t) from equation (13), each curve being for a constant time, t, 


where t changes in increments of 20 hrs beginning at 0.5 hrs. 


Computations are for K = 0.0085 cm} and Q; = 1.75 cal cm hr}, 
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Figure 4c. Ty(Z,t) from equation (14), each curve being for a constant time, t, 


where t changes in increments of 20 hrs beginning at 0.5 hrs. 
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5. (a) Typical salinity (parts per thousand) and temperature (°c) vs 
depth (cm) for winter ice in April. 

Using conditigns for winter ice in (a), computations of the, 
diffusivity a (cm hr ) and thermal conductivity k(cm (CG Tipke 
vs depth (cm) are shown. e 

(c) Typical salinity (parts per thousand) and temperature (C) vs 
depth (cm) for polar ice in April. 
Using conditigns for polar ice sing (Cc); f 
diffusivity a (cm hr ) and thermal conductivity k(cm 
vs depth (cm) are shown. 
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Figure 6. Computations of heat fluxes and ice temperatures for a 12-day period 


following 2300 hrs (LMT) on the 15th of April, 1968, using an 


albedo a = 0.65 and an extinction coefficient K = 0.0175 ae. 
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Figure 7. Computations of heat fluxes and ice temperatures for a 12-day period 


following 2300 hrs (LMT) on the 15th of April, 1968, using 


an albedo a = 0.75 and an extinction coefficient K = 0.034 cm, 
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Figure 8. Computation of heat fluxes and ice temperatures for a 12-day period 


following 2300 hrs (LMT) on the 15th of April, 1968, using an 


albedo a = 0.55 and an extinction coefficient K = 0.0085 pee 
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